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Abstract 

The present paper deals with a next-to-leading order analysis of the 
radiative B — > K*j decay. Working in the PQCD approach developed 
by Szczepaniak et al, we compute the correction, coming from a single 
gluon exchange with the spectator, to the essential form factor. Since 
the branching ratio gets much above the experimental data, although 
in agreement with other theoretical models predictions, we take into 
consideration the effects of complex flavour couplings in the squark 
sector. Finally, we discuss these SUSY implications on the branching 
ratio and direct CP asymmetry values and impose bounds on the 
squark mixing parameter ($23) LR - 
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After the Cabibbo-favoured b — > 57 mode was first reported, in 1993, 
by CLEO 11 [1] and updated in 1995 [2], the exclusive radiative decays, 
B — > K*y and B — > py, as well as the inclusive ones, B — > X s ^y, have be- 
come main targets for both experimental and theoretical investigations. The 
exclusive modes, which are easier to be experimentally investigated [3, 4, 5], 
but less theoretically clear, have been worked out in different approaches. For 
example, the spin symmetry for heavy quarks combined with wave function 
models [6, 7] or the heavy quark effective theory when both b and s are heavy 
[8] have been used. Also, perturbative QCD (PQCD) formalisms, introduced 
for exclusive nonleptonic heavy-to-light transitions, have been extended to 
account for the radiative decays. Recently, detailed analyses of B — > K*y 
and B — > py, in the next-to- leading order (NLO), with the inclusion of hard 
spectator and vertex corrections, have been performed [9-11] and a consis- 
tent treatment, based on a new factorization formula, has been proposed 
[12]. Besides an independent determination of the |Vfd/Vf S | ratio, the b — > 57 
decays are suitable for studying the viability of SUSY extensions of the SM, 
in view of flavour changing neutral currents (FCNC) and CP tests, and for 
imposing constraints on the supersymmetric benchmark scenarios [13, 14]. 

The aim of the present paper is to analyse the B — > K*y decay, in the 
minimal supersymmetric SM (MSSM) context. First, at next-to-leading or- 
der, we compute the hard-spectator correction to the essential form factor. 
In this respect, we employ the PQCD approach developed by Szczepaniak 
et al. [15], for decays dominated by tree diagrams and later extended to 
"penguin" processes [16]. As the branching ratio gets much above the ex- 
perimental data, we make use of the mass insertion method to include, in 
the Wilson coefficients Cr,s, gluino-mediated FCNC contributions. Finally, 
the available data on Br and direct CP asymmetry are used to constrain the 
complex values of the squark mixing parameter (5f 3 ) 



The effective Hamiltonian which describes the B — > K*y radiative decay 
is given by [9, 10] 



where X p = V p bV* s , with p summed over u and c, and Ci, C 7 , C 8 are the 
effective Wilson coefficients at p — m b . The hadronic matrix elements of the 
four-fermion operator and of the electromagnetic and chromagnetic penguin 
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where is the the momentum of the photon and t v is the K* 4-vector po- 
larization. In the heavy quark limit, m b ^> Aq C d, neglecting the corrections 
of order l/m b and a s , one has the following relation among the form factors 
[9] 



m B 



V(0) 



uib + rn K * 



A 1 (0) = 71(0) = T 2 (0)=F x *(0) (5) 



tub + rriK* rriB 

This relation is broken when one includes QCD radiative corrections coming 
from vertex renormalization and hard gluon exchanges with the spectator. 
We recommend [9, 10, 12] for detailed analyses of both factorizable and non- 
factorizable vertex and hard-spectator contributions, involving the operators 
7 , 8 and penguin- type diagrams of 0±. However, it has been stated that 
factorization holds, at large recoil and leading order in l/m& [11] and quan- 
titative tests for proving QCD factorization at the level of power corrections 
have been provided [17]. 
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For a consistent treatment of radiative decays, at next-to-leading order 
in QCD, a novel factorization formula have been proposed in [12]. In this 
approach, the hadronic matrix elements in (1) are written in terms of the 
essential form factor, which describes the long-distance dynamics and is a 
nonperturbative object, and of the hard-scattering kernels, T/ and T/ 7 , in- 
cluding the perturbative short- distance interactions, as 



(6) 



where 77 is the photon polarization. When the dominant contribution comes 
from 7 , we use (4) to write down the decay amplitude as 
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and consequently the branching ratio reads 
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with z = itlk* /%• At next-to- leading order in a s , one has to consider, 
in (6), the contributions to the hard scattering kernels T/ coming from the 
operators 0\ and O s . These have been evaluated in [12] and bring (7) to the 
expression 
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with s c = ml/ml and \i = m&. 

Going further, we add factorizable NLO hard-spectator corrections, to 
the form factor Fk*(0). For a single gluon exchanged with the spectator 
(see Figure 1), we extend the PQCD approach, developed by Szczepaniak 
et al. for heavy-to-light transitions dominated by tree diagrams [15], to the 
so-called penguin processes. 

We evaluate the matrix element of the operator O7, as the following trace 
over spin, flavor and color indices, and integration over momentum fractions 
[16] 
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where Q 2 ~ —(1 — x)(l — y)m 2 B . The B meson wave function 
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contains a strongly peaked distribution amplitude, around a = 
0.072, for X B = 0.38. The K* is described by the wave function 
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(Figure 1) 
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where the light-cone distribution amplitude, tpK*{y), has the following ex- 
pansion in Gegenbauer polynomials [18] 



y^(y) = 6y(l-y)[l + afcf /2) (2y 



l) + afC 3 2 /2 (2y 



1) + ...], (14) 



with Cl ,2 {u) = 3u, Cf /2 (u) = (3/2) (5m 2 - 1), af(m b ) = 0.18 ± 0.05, and 
' c *(rrib) = 0.03 ± 0.03. Performing the calculations in (11) and using the 
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form factors decomposition (4), we identify the spectator contribution to the 
essential form factor as 
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where the if*-mass is neglected. Since we have introduced a cut-off for y — > 1, 
the form factor correction (15) depends on the peaking parameter a and this 
is a main uncertainty in our calculations. For the following input values: 
a s (j2 = Q 2 ) « 0.38, f B = 0.180 GeV, /£. = 0.185 GeV and a = 0.072, we 
get 

F sp (0.072) = 0.1475 (16) 

With the total form factor F K *(0) + F sp (0.072) = 0.38 + 0.1475 in the am- 
plitude (9), the branching ratio gets significantly enhanced to the value 
Br NLO = 6.97 x 10~ 5 . This is comparable to the average theoretical pre- 
diction (7.5 ± 0.3) x 10~ 5 [9, 11, 12], but is much above the experimental 
data: 



Br(B + -> fT+7) 



Br(B° -> K*°-y) = < 



(3.83 ± 0.62 ± 0.22) x 10~ 5 (BaBar [3]) 

(3.76lg;|l ± 0.28) x 10- 5 (CLEO [4]) 

(3.89 ± 0.93 ± 0.41) x 10~ 5 (Belle [5]) 

(4.23 ± 0.40 ± 0.22) x 10~ 5 (BaBar [3]) 

(4.55!8;H ± 0.34) x 10~ 5 (CLEO [4]) 

(4.96 ± 0.67 ± 0.45) x lO" 5 (Belle[5]) 



Moreover, the direct CP asymmetry, defined as 

T(B -> K*j) -T(B -> AT* 7 ) 



a C p = 



T(B -> K*-f) + F(B -> AT* 7 ) 



(17) 



has been predicted by the SM to be acp < 1 0.005 1 , and this disagrees with 
the BaBar and CLEO data, a C p = -0.044 ± 0.076 ± 0.012 (BaBar [3]) and 
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a CP = 0.08 ± 0.13 ± 0.03 (CLEO [4], for the sum of neutral and charged 
B — > K**y decays). 

So, predictions for measurable values of the direct CP asymmetry, in 
agreement with more precise measurements, might have a dominantly new 
physics origin. Following this idea, let us analyse the B — > K*^y decay in the 
MSSM context. Using the mass insertion approximation, [19], we incorpo- 
rate, in the Wilson coefficients C 7 and C 8 , the FCNC SUSY contributions 



where 



4a; 

F (x) = - — -7 1 + 4a; - 5a; 2 + 4a; ln(x) + 2a; 2 ln(x) 
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G (x) = ^ [22 - 20a: - 2a; 2 + 16s ln(x) - a; 2 ln(g) + 9 ln(x)l (19) 
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In (19), a; = m|/m| is expressed in terms of the gluino mass, and 
an average squark mass, rriq. We underline that, in the expressions of 
Cf % (Msusy), we have kept only the left-right squark mixing parameter 
(^23) = (A-bs) / m 1 since, being proportional to the large factor m g /mb, 
will have a significant numerical impact on the branching ratio value. The 
quantities Ab s are the off-diagonal terms in the sfermion mass matrices, con- 
necting the flavours b and s along the sfermion propagators [19]. In these 
assumptions, the total Wilson coefficients, encoding the New Physics, become 

C^ix^] = C 7 (m b ) + C* USY (m b ) , 

Ct tal [x,5] = C 8 (mb) + Ci USY (m b ), (20) 

where C^ SY (m b ) have been evolved from M S usy = m g down to the /x = m h 
scale, using the relations [20] 

Ci USY (m b ) = r)Ci USY (m,), 

Cf SY (m h ) = V 2 Cf USY (m~ g ) + 8 -( V - V 2 )Ci USY (m,) , (21) 
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with 

rj = (a s (mg)/a s (m t )) 2/21 (a s (m t )/a s (m b )) 2/23 ( 22) 

Finally, putting everything together, we replace, in (9), the Wilson coef- 
ficients C 7 and C 8 respectively by Cj otal [x, 5] and Cl° tal [x, 5], the form factor 
F K *(0) by F K *(0) + F sp (a) and, consequently, the branching ratio (8) turns 
into 



Br total = Br SM + SUSY = ^ (1 _ ^2 )3 | F ^ (Q) + p s P {a) 

327T 4 
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One can notice that, for a given x and 5 = pe iip , the total branching ratio is 
depending on three free parameters: a, p, <p, while the direct CP asymmetry 
parameter, (17), is free of the uncertainty a. 

In the next coming discussion, we use the following input parameters: 
m b {m b ) = 4.2 GeV, a = 1/137, \V tb V*\ = 0.0396±0.002, r B o = (1.546 ± 0.018) 
ps, the QCD sum rules analyses result F K *(0) = 0.38, and rriq = 500 GeV. In 
what it concerns the gluino, as its pair production cross section has large can- 
cellations in the e + e~ annihilation, there is hope that the laser-backscattering 
photons will provide a precise gluino mass determination [21]. For a wide 
range of squark masses, a gluino mass of 540 GeV may be measured, with a 
precision of at least ±2 ... 5, at the multi-TeV linear collider at CERN. 



For x taking the values xi = 0.3, x = (540/500) 2 and x g = 3 (where 1(g) 
comes from rrig less (greater) than rriq), and imposing the BaBar constraint 

[3] 

-0.17 < a CP < 0.082, (24) 

we draw, in Figure 2, the contour plots of constant Br total (the dashed lines) 
and a C p (the solid lines). When {p, ip} £ [0, 0.03] x [— n/2, vr/2], we get, for 
the world average branching ratio data, over the B ± and B° decay modes, 
Br exp (B ± -> K M± i) = (4.22 ± 0.28) x 10" 5 , three dashed lines, with in- 
creasing thickness, as x goes from xi to x g . Correspondingly, for acp, we get 
three pairs of solid curves: the lower ones, for acp = —0.17, and the upper 
ones, for acp = 0.082. These solid contours close inside the values of direct 
CP asymmetry which do not agree with (24). Now, we are able to put strong 



(Figure2) 
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constraints on ($23) LR i by looking at the segments of the .Br-plots outside 
the solid contours, for each x. We notice that, for rrig > m,g, all the negative 
phases, with suitable p's, can accommodate both the relation (24) and the 
branching ratio data. 

In Figure 3, we represent, the Br total (in units of 10~ 5 ) and a>cp (in 
units of 0.1), with respectively dashed and solid lines, as functions of <p, for 
x = 3. As p takes the following values: p E {0.005, 0.01, 0.015, 0.02}, we get 
four pairs of curves, with increasing thickness. The horizontal dashed line 
corresponds to the average Br data, while the horizontal solid ones stand for 
the constraint (24). 

Finally, let us perform a numerical analyses, for x = xq, and increas- 
ing p, starting with p = 0.005. As ip E [-8tt/16, -4tt/16] U [3tt/16, 7tt/16], 
the Br total and the direct CP asymmetry are inside the ranges 10 5 x Br total E 
[8.3, 3.3] and [3.1, 8.1] and, respectively, a C p E [-0.054, -0.093]U[0.096, 0.062], 
accommodating data and other theoretical models predictions. When p goes 
to bigger values, the two <p ranges, constrained by the allowed branching ra- 
tios, get closer and acp moves toward much bigger values. For example, for 
p = 0.01 and (p E [-6tt/16, -47r/16]U[37r/16, 5tt/16], one gets 10 5 x Br total E 
[8.1, 3.6] and [3.2, 7.62] and, respectively, a CP E [-0.1, -0.16] U [0.21, 0.12] 
and we notice that only the negative (^-values lead to acp inside the BaBar 
constraint (24). For p = 0.015 and tp E [-4tt/16, -2tt/16] U [tt/16, 3tt/16], 
the values 10 5 x Br total E [7.8, 3.6] and [3.3, 7.3] are compatible with a mea- 
surable acp ~ ±0.12. Starting with p = 0.02, the predictions for branching 
ratio are above data and other theoretical estimations, the minimum value 
being Br total = 9 x 10~ 5 , for ip = 0, while the corresponding asymmetry is 
a CP = 0.003. 

In the present paper, we have analysed the radiative B — > K*^ de- 
cay, in a combined PQCD and SUSY framework. First, we have used the 
PQCD approach, developed by Szczepaniak et al. [15] and extended to "pen- 
guin" processes [16], to compute the hard-spectator contribution, F sp (a), 
to the essential form factor F K *(0). For the peaking parameter in the B 
wave function a = 0.072 and F K *(0) = 0.038, the branching ratio becomes 
Br NLO = 6.97 x 10~ 5 , which is above the experimental data, while the di- 
rect CP asymmetry predicted by the SM lies much below [3, 4]. In order to 



(Figure 3) 
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find an agreement, we extend our analyses by including, in the Wilson coef- 
ficients CV,8, the SUSY contributions coming from squark mixing parameter 



a, on three (SUSY) parameters: x,p,<p, while acp is free of the uncertainty 
coming from the form factors. Using the graphs displayed in Figures 2 and 
3, one is able to find out allowed ranges for the mass insertion parameter 
($23) LR - As an example, for x = (540/500) 2 , the world average branching 

ratio, Br exp = 4.22 x 10~ 5 , can be accommodated for {p, ip} = |o.005, — || j 

or jo. 01, — |||. The corresponding asymmetries, acp = —0.085 and respec- 
tively acp = —0.147, are inside the BaBar constraint (24). 
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FIGURE CAPTIONS 



Fig.l. The Feynman contributing diagrams in the hard scattering amplitude 
T M . The gluon and photon are respectively represented by dotted and dashed 
lines. 

Fig. 2. Contour plots of total branching ratio fitting the world average data, 
in units of 10~ 5 , (the dashed lines) and the BaBar constraint (24) on direct 
CP asymmetry (the solid lines), as functions of p and (p. The thickness of 
contours is increasing as x takes the values: x = {0.3, 1.16, 3}. The solid 
curves close inside the values of acp which disagree with the constraint (24). 
Fig. 3. The total branching ratio, in units of 10~ 5 , (the dashed lines) and 
10 x a CP (the solid lines), as functions of </?, for x = 3. The thickness 
of plots increases as p is respectively: p = {0.005, 0.01, 0.015, 0.02}. The 
horizontal dashed line corresponds to the world average branching ratio and 
the horizontal solid lines are for the constraint (24). 
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